Evolution of electronic structure of doped Mott insulators 
- reconstruction of poles and zeros of Green's function 
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We study evolution of metals from Mott insulators in the carrier-doped 2D Hubbard model using a cluster 
extension of the dynamical mean-field theory. While the conventional metal is simply characterized by the Fermi 
surface (pole of the Green function G), interference of the zero surfaces of G with the pole surfaces becomes 
crucial in the doped Mott insulators. Mutually interfering pole and zero surfaces are dramatically transferred 
over the Mott gap, when lightly doped holes synergetically loosen the doublon-holon binding. The heart of the 
Mott physics such as the pseudogap, hole pockets, Fermi arcs, in-gap states, and Lifshitz transitions appears as 
natural consequences of this global interference in the frequency space. 

PACS numbers: 71.10.Hf; 71.30.+h; 74.72.-h 



Electronic structures of doped Mott insulators have inten- 
sively been debated since the discovery of high-T c cuprates. 
While overdoped metals behave like a conventional Fermi liq- 
uid, an outstanding issue is how such metals evolve into the 
Mott insulator. Experimental studies on the doped cuprates 
revealed anomalous metallic behaviors in the proximity to the 
Mott insulator. Among them angle-resolved photoemission 
spectroscopy (ARPES) found arclike spectra around the nodal 
points of the momenta k = (±§, ±§), with a pseudogap in 
the antinodal region of the Brillouin zone [ 1]. Moreover the 
high-resolution ARPES [2] and quantum oscillations of resis- 
tivity observed under magnetic fields [3] suggested the exis- 
tence of a small closed Fermi surface, called the pocket, in 
contrast to the large surface expected in the normal Fermi liq- 
uid. These imply emergence of non Fermi liquids under a 
radical evolution of electronic structure upon doping. 

Diverse theoretical proposals were made for the Mott 
physics of the 2D Hubbard model: The pseudogap was re- 
produced in various theoretical frameworks 6]. The ex- 
act diagonalization (ED) studies [7] further suggested a spec- 
tral weight transfer with doping from the upper Hubbard band 
(UHB) to the upper edge of the lower Hubbard band (LHB), 
which created in-gap states. The arc structure of the spectral 
weight was reproduced in the cluster perturbation theory |8|], 
the dynamical cluster approximation [|9|,LL0|], and the cellular 



dynamical mean-field theory (CDMFT) M& lilt LL2l LLJ • Hole- 
pocket Fermi surfaces around the nodal points were suggested 
in a phenomenology [14], the CDMFT 11211 . and a variational- 
cluster approach 111511 . A topological change of the Fermi sur- 
face, i.e., Lifshitz transition [16], from electronlike to holelike 
surfaces due to a correlation effect was analyzed 
Lifshitz transition to electron pockets was also found upon 
electron doping [19]. In spite of these achievements, a coher- 
ent picture of the Mott physics has not emerged yet. 

In general, poles of the single-particle Green function 
G(k, uj) are known to define the Fermi surface at the fre- 
quency uj = 0. While Re G changes its sign across a pole 
through ±oo, Re G may also change the sign across a zero 
defined by G = 0. Recent studies suggest, in addition to the 



insulator, the self-energy E(k, uj) inside the gap diverges on a 
specific surface in the k-uj space, which defines a zero surface 
of G. Since the Fermi surface disappears in the Mott insula- 
tor while the zero surface appears instead, the evolution of the 
zeros makes crucial effects at low doping and is imperative in 
understanding the Mott physics. 

In this letter we study the 2D Hubbard model using the 
CDMFT+ED method 11211 . and clarify the doping evolution of 
poles and zeros of G in the \t-uj space. By figuring out recon- 
struction of interfering poles and zeros in a wide uj range, frag- 
mentary features are integrated into a coherent understanding 
of the Mott physics at and around uj = 0, namely, the hole- 
pocket Fermi surface, Fermi arc, pseudogap, in-gap states, and 
Lifshitz transitions. Starting from a high-energy structure in 
the scale of the Mott gap, we zoom in lower-energy hierar- 
chy, in a clear and visible fashion, revealing two structures 
induced by hole doping. One is in-gap states formed above 
the occupied states separated by a pseudogap. The other is 
hole pockets in the nodal directions, which appear at the top 
of the occupied states below the pseudogap. The pseudogap 
is described by a zero surface crossing the Fermi level. The 
density-of- states analysis suggests that these structures result 
from a dramatic transfer of the weight over the Mott gap: This 
is caused by tiny doping, which loosens the doublon-holon 
binding. Then the zero surface interferes with the poles and 
pushes down the pole surface near the zeros below the Fermi 
level, leaving hole pockets around (±f , ±§), which are more 
stabilized by the next-nearest transfer. With further doping, 
the Fermi surfaces undergo at least two continuous Lifshitz 
transitions before reaching a normal Fermi liquid with an elec- 
tronlike Fermi surface. 

The Hubbard Hamiltonian on a square lattice reads 



H = Ye(k)cl a c^ 



(1) 



where c\^ a (c\. a ) destroys (creates) an electron of spin a with 
momentum k, and is a spin density operator at site i. U 
represents the onsite Coulomb repulsion, ji the chemical po- 
tential, and e(k) = — 2t(cos k x ± cos k y ) — At' cos k x cos k y , 



poles, zeros of G play important roles iSEllll. In the Mott where t (t f ) is the (next-)nearest-neighbor transfer integral 
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FIG. 1 : The k-u; structure of poles and zeros of G for t' — 0, U = 8t 
at (a) £=0 and (b) 0.09. Poles (zeros) are plotted in green (red). Oc- 
cupied regions for an electron are filled with aqua. The transparency 
of pole (zero) surfaces reflects the imaginary part of the (inverse) 
self-energy (more transparent for larger imaginary part), (c) and (d) 
show the density of states per spin at 5 = and 0.09, respectively. 

In the CDMFT, the infinite lattice of the model CQ) is self- 
consistently mapped onto an A^-site cluster embedded in in- 
finite bath sites, which define an N c x N c dynamical mean- 
field matrix go(iw n ) for the cluster in terms of Matsubara fre- 
quency uo n . We employ the ED method with the Lanczos algo- 
rithm to solve the cluster problem at zero temperature, where 
go is fitted with a finite number of bath parameters and a ficti- 
tious temperature 1/(3 is introduced. We adopt 2 by 2 cluster 
(N c = 4) coupled with 8 bath sites and (3 = 100/t unless oth- 
erwise mentioned. When a self-consistency loop converges, 
we obtain an N c x N c self-energy matrix E c (icc; n ), whose el- 
ements are defined on sites within the cluster. Then we inter- 
polate cluster quantities in the momentum space, to obtain the 
original infinite-lattice ones. To study poles and zeros of G si- 
multaneously, after careful examination, we employ a suitable 
scheme of the cumulant periodization 11211 , where the clus- 
ter cumulant M c = (ico n + fi — E c ) _1 is interpolated with a 
Fourier transformation cut by a cluster size N c , 

1 Nc 

M(k,^ n ) = — ^ [McWliie*'^-^). (2) 
c i,i=l 

After M is obtained, G and E on the original lattice are cal- 




FIG. 2: (a)(b) Close-ups of Fig. Q2b) at two different lower-energy 
scales. Blue grade levels in (b) depict rigid-band shifts of the chem- 
ical potential. (c)(d)(e) r(k) = Re[G f_1 (k, 0)] in the first quadrant 
of the Brillouin zone for t' = and U = St at S =0.09, 0.11 and 
0.14, respectively, r = (=Loo) defines the Fermi surface (the zero 
surface of G). We set 77 = 10 _4 £ for these data. 



culated with G = [M _1 - e(k)]" 1 and E = iu n + /1 - M ~\ 
respectively. The Lanczos-ED method allows us to calculate 
G at real frequencies directly through a continued-fraction ex- 
pansion 112 ill , where a Matsubara frequency iuo n is replaced by 
u + if] with a small positive factor 77. Although we should be 
careful on the limitation of the present size of the cluster, this 
combination of the methods is appropriate as the state of the 
art in exploiting zero and pole structures on equal footing in a 
wide spectrum range. 

We start from electronic structure of the Mott insulator at 
half filling (5 = 0; 5 is the hole doping rate from half filling) 
for t' = and U/t = 8. While an insulating phase caused by 
some symmetry breaking may show a similar feature, to cap- 
ture the essence of the Mott physics, we assume a paramag- 
netic solution in the CDMFT. Figure [TJa) shows the structure 
of poles (green surfaces) and zeros (red surfaces) of G in the 
k-cj space. Here uo is measured from the upper edge of the 
occupied states B22h to compare with the doped case below. 
The zero surface around < uj < 3t (the most distinct red 
surface) represents the Mott gap, which separates UHB and 
LHB displayed by entangled pole and zero surfaces. The zero 
surface at the Mott gap touches the pole surfaces of the UHB 
(LHB) around (0,0) [(tt, tt)] at u ~ 3t (u ~ 0). The distinct 
pole surface in the UHB persists up to co ~ 4t, and another 
distinct one appears around bt < uj < It. In accordance with 
these pole surfaces, there are large weights in the density of 
states (DOS) per spin [Fig.QJc)]. In other area, poles are dis- 
rupted due to a large imaginary part of the self-energy, making 
incoherent contributions. 

We next dope holes. Comparing the result at 5 = 0.09 in 
Fig- Gib) with Fig. [2a), we see that high-energy > 5f) 
structure does not change appreciably while the lower-energy 
(| a; I < bt) structure does: The pole surface at the bottom of 
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FIG. 3: (Color online) (a) DOS at S = 0.09 [the same as Fig. Ed)] 
and the definitions of weights Wi and Wi. (b) S dependence of Wi 
and W2 compared with the weights estimated from a simple analysis 
given by dashed lines (see the text). We set 77 = 10 _2 £ for these data. 



UHB just above the Mott gap ~ 3t in Fig. [Ha) is transferred 
to the lowest-energy pole surface in the piles above uj ~ 0.2t, 
squashing the zero surface just below it, as in Fig.QJb) and as 
will be seen in Figs. [2a) and (b). In addition, one of vague 
zero surfaces in the UHB piles in Fig. Oa) turns to a distinct 
one at around co = 4t in Fig. \Hb). Corresponding to these 
reconstructions, the DOS at co ~ 3t in Fig. [He) is transferred 
to energies lower than co ~ 2t as shown in Fig.QId). 

Let us focus on a low-energy structure of poles and zeros. 
Figures [3a) and (b) provide enlarged views of Fig. [Tib) near 
the Fermi level, in two different energy scales. A prominent 
feature is a small gap, described by a zero surface cutting the 
Fermi level. The gap has a small width of ~ 0.2t, and is dis- 
tinguished from the larger gap between the doped hole states 
and the UHB as shown in DOS in Fig. Htd). A comparable 
pseudogap was found in earlier CDMFT studies Il6lll2ll. 

As the zero surface around (tt, tt) extends to the region 
uo < 0, near the zero surface, the pole surface is pushed down 
below the Fermi level [Fig. Ob)], because a zero and a pole 
surfaces cannot cross each other. Meanwhile in the regions 
far away from zeros, the energy of poles increases with |k|, 
reflecting the original dispersion e(k). Hence, along the di- 
rection from (0,0) to (7r, tt) for example, the pole energy ini- 
tially increases, crossing the Fermi level, and then turns down 
around (f , §), crossing the Fermi level again. As a conse- 
quence, a hole pocket is formed around (f , f ), accompany- 
ing a zero surface around (tt, tt), as found previously in Ref. 
12l and also shown in Fig.O c )- 

The metal-insulator transition occurs when the hole pock- 
ets shrink and vanish, indicating a topological nature of this 
transition [23]. Since the quick evolution for S < 0.01 eluci- 
dated below suggests a proximity to the first-order transition, 
the marginal quantum criticality J23ll may have relevance. 

To get further insight into the underlying physics of the dra- 
matic restructuring, we study the doping evolution of DOS 
systematically. We calculate the weight integrated over en- 
ergy ranges of < uo < 2. St (in-gap weight), Wi, and that 
over 2. St < co < 00 (UHB weight), W2 [see Fig. Oa)]. Figure 
Ob) depicts S dependences of W\ and W2, where W\ (W2) 
monotonically increases (decreases) with S. 



As shown in Fig. Ob), the in-gap weight W\ shows good 
agreement with S + rid, where rid = (^it n *l) i s me doublon 
density. It is natural to have the weight proportional to S in the 
rigid band picture because it comes from the holes doped into 
the LHB, but there is an additional contribution scaled by rid. 
This is interpreted as a strong-correlation effect as follows. 
Neglecting dynamical fluctuations, we can roughly estimate 
the single-occupancy weight as Wo = ^^2 a (^a — ^d) = 
|(1 — S — 2nd), where n G = (n^) (a =j, j) is the spin 
density. Now, one electron added to these singly occupied 
states makes a "doublon" and contributes to the DOS of the 
UHB. This contribution should agree with W2 . Since the DOS 
below the Fermi level has a weight |(1 — 5), the remaining 
weight, which should be equal to Wi, becomes 1 — ^(1 — 5 — 
2rid) — ^(1 — S) = S + rid. This simple analysis explains well 
the overall behavior of Wi and W2, as seen in Fig. Ob). 

A nontrivial finding, however, is that the weight rid in Wi 
is very quickly transferred from the UHB to the top of LHB 
by tiny doping S < 0.01 to the Mott insulator. Actually, the 
above simple analysis breaks down at very small S, where W\ 
eventually and obviously vanishes toward S = as shown in 
Fig. Ob). At S = 0, this weight rid is included in the UHB 
weight W2 and is interpreted as an electron excitation added 
to a holon tightly bound with a doublon. The number of dou- 
blons and holons should be the same in the Mott insulator and 
they are tightly bound with the binding energy equal to the 
Mott gap ~ 3t. Therefore one electron added to a hole site re- 
quires dissolution of the bound state and overcoming the Mott 
gap. The quick transfer of rid from W2 to W\ means that the 
doublon-holon binding energy is quickly reduced from ~ 3t 
to ~ 0.2t presumably because of an efficient screening of the 
doublon-holon interaction by the doped holes. Even with this 
screening, the binding energy though small still survives as 
the pseudogap. This quick reduction is promoted by a posi- 
tive feedback, because the dissolved bound pairs further join 
in screening other bound pairs. This constitutes a mechanism 
for the drastic restructuring in the Mott physics. 

The next issue is how the hole-pocket Fermi surface ac- 
companied by zeros of G in Fig. Oc) evolves into a single 
electronlike Fermi surface expected in the normal Fermi liq- 
uid as S increases further. Figures Oc), (d), and (e) depict 
r(k) = RelG-^k,*))] at S = 0.09, 0.11, and 0.14, respec- 
tively, showing how Fermi (zero) surfaces at r = (±00) 
evolve. As 5 increases from 0.09, the hole pocket continues to 
expand until it touches the Brillouin zone boundary (\k x \ = ir 
or I fey I = 7r). When it touches the boundary, it changes into 
two concentric Fermi surfaces around (tt, tt) through a Lifshitz 
transition, while the zero surface remains around (tt, tt) [Fig. 
Od)] . Further doping enlarges the unoccupied region (blue re- 
gion) sandwiched by the two concentric Fermi surfaces. Then, 
the smaller Fermi surface around (tt, tt) merges with the zero 
surface and they are annihilated in pair, leaving a large hole- 
like Fermi surface, which almost simultaneously transforms 
into a normal electronlike one through another Lifshitz transi- 
tion [Fig.Oe)]. For S > 0.14 only a large electronlike Fermi 
surface is found. This evolution of Fermi and zero surfaces is 
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FIG. 4: (Color online) A(k,0) = -Im[C(k,cj + irj)]/ic for t' 
-0.2t, U = 12t, and 6 = 0.09 with (a) 77 = 10" 4 t and (b) 77 
0.05£. Here we set the fictitious temperature 1//3 = t/200. 



understood in a rigid band picture for the electronic structure 
in Fig. Ob), as drawn with blue grade levels, if we assume that 
the hole doping only lowers the Fermi level without changing 
the structure of poles and zeros. Detailed inspection of the two 
Lifshitz transitions shows continuities of the electron density 
as a function of fi and hence continuous transitions. 

Apparently, the Luttinger sum rule B24I1 is violated for S < 
0.11 while roughly satisfied for S > 0.14. Since the Luttinger 
theorem assumes an adiabatic continuity, there is no reason 
that the rule should be satisfied beyond the Lifshitz transition 
around S = 0.14. Other numerical calculations also pointed 
out the violation at small S ifT^L [TtL [Tgh ■ 

Next we consider effects of t' . Results for 5 = 0.09 (not 
shown) exhibits that t' < elevates (lowers) the energy of 
poles around the nodal direction [around (tt, 0) and (0, 7r)], as 
expected from the original dispersion e(k). The enhancement 
around the nodal direction more stabilizes the hole-pocket 
structure while the reduction around (tt,0) and (0, 7r) more 
stabilizes a single holelike Fermi surface. Quantitative modi- 
fications caused by nonzero t' will be reported elsewhere. 

Finally we discuss how a Fermi arc emerges in the lightly 
doped region. Stanescu and Kotliar found that the zero sur- 
face reduces the spectral weight from the neighboring Fermi 
surface, resulting in a Fermi arc [12]. In the present Lanczos- 
ED method, however, one has to take the limit of 77 — > + in 
principle, which leads to ^-function peaks in ImE at the zero 
surfaces of G. In this limit there is no contribution from the 
zeros to the Fermi surface region, namely, the closed hole- 
pocket structure remains robust in the spectral weight and the 
Fermi arc does not show up, as demonstrated in Fig. Ua) for 
a sufficiently small 77. In real situations of experiments, how- 
ever, ImE is broadened through various extrinsic factors, such 
as temperature, impurity scatterings and phonons W25I1 . These 
effects are phenomenologically represented by a broadening 
to ImH through a nonzero 77. In fact, for larger 77, the singu- 
larities of ImE are smeared out and the spectral weight close 
to the zeros are suppressed as shown in Fig. EJb). Thus the 
Fermi arc is reproduced by introducing a phenomenological 
broadening factor 77 in the present result. 

To summarize, we have shown that key elements of the 
Mott physics such as in-gap states, hole pocket, Fermi arc, 
pseudogap, and Lifshitz transitions result from global recon- 
structions and interferences of pole and zero surfaces with 



doping to the Mott insulator. The reconstructions are caused 
by a drastic relaxation of the doublon-holon binding at tiny 
doping. The criticality at 5 — > and a mechanism of the re- 
laxation, though we have given a qualitative picture, should 
further be clarified in future studies. It is also desired to con- 
firm our results in calculations for larger clusters. 
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